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The electron spin is a very promising candidate for a 
solid-state qubit [1]. Major experimental breakthroughs 
have been achieved in recent years using quantum dots 
formed in semiconductor heterostructures based on GaAs 
technology d, IE fill]- In such devices, the major sources 
of spin decoherence have been identified as the spin-orbit 
interaction, coupling the spin to lattice vibrations 10,0, Hi, 
and the hyperfine interaction of the electron spin with the 
surrounding nuclear spins 02, Ell [U Gjl [3 Qj]. There- 
fore, it is desirable to form qubits in quantum dots based 
on other materials, where spin-orbit coupling and hyper- 
fine interaction are considerably weaker [15]. It is well 
known that carbon-based materials such as nanotubes or 
graphene are excellent candidates. This is so because spin- 
orbit coupling is weak in carbon due to its relatively low 
atomic weight, and because natural carbon consists pre- 
dominantly of the zero-spin isotope 12 C, for which the hy- 
perfine interaction is absent. Here we show how to form 
spin qubits in graphene. A crucial requirement to achieve 
this goal is to find quantum dot states where the usual val- 
ley degeneracy is lifted. We show that this problem can be 
avoided in quantum dots with so-called armchair bound- 
aries. We furthermore show that spin qubits in graphene 
can not only be coupled (via Heisenberg exchange) be- 
tween nearest neighbor quantum dots but also over long 
distances. This remarkable feature is a direct consequence 
of the Klein paradox being a distinct property of the quasi- 
relativistic spectrum of graphene. Therefore, the proposed 
system is ideal for fault-tolerant quantum computation, 
and thus for scalability, since it offers a low error rate 
due to weak decoherence, in combination with a high error 
threshold due to the possibility of long-range coupling. 

Only very recently, the fabrication of a single layer of 
graphene and the measurement of its electric transport proper- 
ties have been achieved lfl6l ITtI [l8tl . Two fundamental prob- 
lems need to be overcome before graphene can be used to 
form spin qubits and to operate one or two of them as pro- 
posed in Refs. llllSt]: (i) It is difficult to create a tunable quan- 
tum dot in graphene because of the absence of a gap in the 
spectrum. The phenomenon of Klein tunnelling makes it hard 
to confine particles fl9l I20I l2lTl . (i i) Due to the valley de- 
generacy that exists in graphene l22l[23il24ll . it is non-trivial 
to form two-qubit gates using Heisenberg exchange coupling 
for spins in tunnel-coupled dots. Several attempts have been 
made to solve the problem (i) such as to use suitable transverse 
states in graphene ribbons to confine electrons 112511 . to com- 
bine single and bilayer regions of graphene l26ll . or to achieve 
confinement by using inhomogeneous magnetic fields l27ll . 
The problem (ii) has not been recognized up to now. Here 
we propose a setup which solves both problems (i) and (ii) 
at once. Similar to Ref. ll25ll we choose to confine electrons 



by using suitable transverse states in a ribbon of graphene, 
cf. Fig. [T] In particular, we assume semiconducting arm- 
chair boundary conditions to exist on two opposite edges of 
the sample. It is by now feasible to experimentally identify 
ribbons of graphene with specific boundaries on the atomic 
scale. These are preferably of zigzag or of armchair type. For 
an experimental realization of our proposal, one would have to 
look for ribbons with semiconducting armchair boundaries. It 
is known that in such a device the valley degeneracy is lifted 
||28[ l29ll . which is the essential prerequisite for the appear- 
ance of Heisenberg exchange coupling for spins in tunnel- 
coupled quantum dots (see below), and thus for the use of 
graphene dots for spin qubits. We show below that spin qubits 
in graphene can not only be coupled between nearest neighbor 
quantum dots but also over long distances. This long-distance 
coupling mechanism makes use of conduction band to valence 
band tunnelling processes and is, therefore, directly based on 
the Klein paradox in graphene 1201 12111 . 

We now discuss bound-state solutions in our setup, which 
are required for a localized qubit. We first concentrate on a 
single quantum dot which is assumed to be rectangular with 
width W and length L, see Fig.[TJ The basic idea of forming 
the dot is to take a ribbon of graphene with semiconducting 
armchair boundary conditions in ^-direction and to electri- 
cally confine particles in y-direction. The low energy prop- 
erties of electrons (with energy s with respect to the Dirac 
point) in such a setup are described by the 4x4 Dirac equation 
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¥ + eViyfV = sV, (1) 



where the electric gate potential is assumed to vary stepwise, 
V(y) = Vgate in the dot region (where < y < L), and 
V(y) = Vbarrier in the barrier region (where y < or y > L). In 
Eq. (Q~|i, cr x and <r y are Pauli matrices, h is Planck's constant 
divided by 2jt, v the Fermi velocity, and e the charge of an 
electron. The four component spinor envelope wave function 
¥ = Q¥ ( f>,^f > ,- X ¥f' ) ,-Wf" > ) varies on scales large com- 
pared to the lattice spacing. At this point, we are only inter- 
ested in the orbital structure of the wave function. The spin 
degree of freedom is neglected until the final part, where we 
discuss the Heisenberg exchange coupling for spins in tunnel- 
coupled quantum dots. In the wave function *P, A and B refer 
to the two sublattices in the two-dimensional honeycomb lat- 
tice of carbon atoms, whereas K and K' refer to the vectors 
K and K' in reciprocal space corresponding to the two valleys 
in the bandstructure of graphene. The appropriate semicon- 
ducting armchair boundary conditions for such a wave func- 
tion have been formulated in Ref. B28I1 and can be written as 
(a = A,B) 



(2) 
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FIG. 1 : Schematic of a graphene double quantum dot. Each dot 
is assumed to have length L and width W. The structure is based 
on a ribbon of graphene (grey) with semiconducting armchair edges 
(white). Confinement is achieved by tuning the voltages applied to 
the "barrier" gates (blue) to appropriate values such that bound states 
exist. Additional gates (red) allow to shift the energy levels of the 
dots. Virtual hopping of electrons through barrier 2 (thickness d) 
gives rise to a tunable exchange coupling J between two electron 
spins localized in the left and the right dot. The exchange coupling 
is then used to generate universal two-qubit gates. 



corresponding to a width W of the ribbon shown in Fig. 
where W is not an integer multiple of three unit cells. The + 
signs in Eq. (0 (as well as in Eq. (0 below) correspond to the 
two possible choices of a number of unit cells that is not an 
integer multiple of three. The full set of plane wave solutions 
of Eq. (0 is readily determined ll29ll . It is well known that 
the boundary condition (0 yields the following quantization 
conditions for the wave vector k x = q„ in x-direction ll28l l29tl 



q n = (n± l/3)n/W, n e 



(3) 



An explicit form of the corresponding wave functions is pre- 
sented in App. A and App. B. The level spacing of the modes 
(0 can be estimated as Ae w Hv7t/3W, which gives Ae as 
30 meV, where we used that v « 10 6 m/s and assumed a quan- 
tum dot width of about W ~ 30 nm. Note that Eq. (0 also 
determines the energy gap for excitations as E gap = 2hvqo- 
Therefore, this gap is of the order of 60 meV, which is un- 
usually small for semiconductors. This is a unique feature 
of graphene that will allow for long-distance coupling of spin 
qubits as will be discussed below. 

We now present in more detail the ground-state solutions, 
i.e. n = in Eq. (0. The corresponding ground-state energy 
e can be expressed relative to the potential barrier V = Vbamer 
in the regions y < and y > L as e = e Vbamer ± ^ v (?n + ^ 2 )^ 2 - 
Here, the ± sign refers to a conduction band (+) and a valence 
band (-) solution to Eq. ([]]). For bound states to exist and to 
decay at y — > +oo, we require that fivqo > \e - e Vbamer I > which 
implies that the wave vector k y = k in ^-direction, given by 



is purely imaginary. In the dot region (0 < y < L), the 
wave vector k in y-direction is replaced by k, satisfying e = 
eVgate ± h v (lo + k 2 ) 1 ^ 2 - Again the ± sign refers to conduction 
and valence band solutions. In the following, we focus on 
conduction band solutions to the problem. 

Since the Dirac equation (03 implies the continuity of the 
wave function, the matching condition at y — and y — L 
allows us to derive the transcendental equation for e 
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with zo,k = (<7o + ft)/(9o + k 2 ) 112 . Eq. (0 determines the al- 
lowed energies e for bound states. In order to analyze the 
solutions to Eq. {3), we distinguish two cases, one where k 
is real, and the other, where k is purely imaginary. The two 
cases are distinguished by the condition \e - <?V ga tel > hvqo 
and |e - eVg at e| < tivqQ, respectively. Furthermore, we assume 
that Vgate + Vbamer, i.e., zo,k + z j- If we relax this assump- 
tion, we can show that for the case zo,k = Zoi on ty a single 
solution to Eq. (0 exists, namely Zqi = 1, which implies that 
k — 0. The corresponding wave function to this solution van- 
ishes identically (see App. A for further details). In the case 
where k is purely imaginary, there is no bound-state solution. 
This is due to the fact that such a solution would have to exist 
directly in the bandgap. We now analyze solutions for real k. 
In the corresponding energy window 



\S - eVgatel > hvqo >\S~ eVbarrierl, 

we can simplify Eq. (0 considerably, obtaining 



tan(fcL) = 



tivk yj '(fivqa) 2 - (s ■ 



€ Vbarrier)*^ 



(£ - eVbarrierXs - e Vg ate ) - (frvq ) 2 



(6) 



(7) 



k = iyjql- ((s- eVba 
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We show a set of solutions to Eq. (0 for a relatively short 
dot (qoL = 2) as well as a longer dot (qoL = 5) in Fig. 
The number of bound states (for n — 0) is maximal if 
AV = Vbamer - V ga te is exactly as large as the size of the gap 
£gap = 2Hvqo, then A^ max = |" y/&qoL/7r\, where \x\ is the in- 
teger just larger than x. The level spacing associated with the 
allowed solutions of Eq. (0 increases as L decreases and has 
a rather complicated parameter dependence. It can, however, 
be estimated to be of the order of Ae * /iv7r/max{W, L), which 
is in the energy range of a few tens of meV as mentioned be- 
low Eq. (0. In Fig. we show the energy bands of a single 
dot and two neighboring barrier regions as well as a double 
dot setup with three barrier regions. The double dot case il- 
lustrates how we make use of the Klein paradox to couple two 
dots. 

A particular example of a wave function is shown in 
Fig.0 It is a ground-state solution under the parameter choice 
e(Vbamei- - Vg a te) = 0.6/zvgo, and qoL = 2 (indicated by the ar- 
row in Fig. 0. The weight of the wave function on the A 
and B lattice sites is different, however, the integrated weight 
is the same as required by the normalization condition f28ll . 
Ground-state solutions (i.e. the lowest lying (red) lines in 
Fig. have no nodes in the dot region - similar to the cor- 
responding problem of confined electrons that obey the non- 
relativistic Schrodinger equation. Excited-state solutions in 
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FIG. 2: Bound state solutions for two different dot sizes. Bound- 
state solutions of a relatively long (qoL = 5, left panel) and a shorter 
(qoL = 2, right panel) quantum dot are shown. The diagonal straight 
lines mark the area in which bound-state solutions can occur. The 
arrow marks the solution for which the wave function is plotted in 
Fig.H 



parameter regions in which they exist, do have nodes in the 
dot region, which is shown in Fig. 8 of App. C. 

We now turn to the case of two coupled graphene quantum 
dots, separated by a potential barrier, as sketched in Fig. Q] 
each dot filled with a single electron. It is interesting to ask 
whether the spins S, of these two electrons (i - 1,2) are 
coupled through an exchange coupling, H exc h — JSi ■ S2, in 
the same way as for regular semiconductor quantum dots [9], 
because this coupling is, in combination with single-spin ro- 
tations, sufficient to generate all quantum gates required for 
universal quantum computation [lfl. The exchange coupling 
is based on the Pauli exclusion principle which allows for 
electron hopping between the dots in the spin singlet state 
(with opposite spins) of two electrons, but not in a spin triplet 
(with parallel spins), thus leading to a singlet-triplet split- 
ting (exchange energy) J. However, a singlet-triplet split- 
ting J + only occurs if the triplet state with two elec- 
trons on the same dot in the ground state is forbidden, i.e., 
in the case of a single non-degenerate orbital level. This is 
a non-trivial requirement in a graphene structure, as in bulk 
graphene, there is a two-fold orbital ("valley") degeneracy of 
states around the points K and K' in the first Brillouin zone. 
This valley degeneracy is lifted in our case of a ribbon with 
semiconducting armchair edges, and the ground-state solu- 
tions determined by Eq. (0 are in fact non-degenerate. The 
magnitude of the exchange coupling within a Hund-Mulliken 
model is U J = (-U H + (U 2 H + 16?^) 1/2 )/2 + V, where t 
is the tunnelling (hopping) matrix element between the left 
and right dot, U is the on-site Coulomb energy, and V is the 
direct exchange from the long-range (inter-dot) Coulomb in- 
teraction. The symbols f# and Uh indicate that these quan- 
tities are renormalized from the bare values t and U by the 
inter-dot Coulomb interaction. For t <K U and neglecting 
the long-ranged Coulomb part, this simplifies to the Hub- 
bard model result J = 4t 2 /U where t is the tunnelling (hop- 
ping) matrix element between the left and right dot and U 
is the on-site Coulomb energy. In the regime of weak tun- 
nelling, we can estimate t « e J x ¥^ L (x,yy\ l R(x,y)dx dy, where 
y) - y ± id + L)/2) are the ground-state spinor 




FIG. 3: Energy bands for single and double dot case, (a) Energy 
bands for two barrier regions and a single dot. The red area marks a 
continuum of states in the valence bands and the blue area marks a 
continuum of states in the conduction bands. In the barrier regions, 
we indicate the energy bands of the quantized modes due to trans- 
verse confinement. All modes are non-degenerate solutions in valley 
space. They come pairwise in a sense that always two of them are 
separated by a distance hvq in energy space. In the figure, this is il- 
lustrated for the energy levels corresponding to wave vectors qo and 
q_i as well as q t and q_ 2 - In the dot region, the electric confinement 
in longitudinal direction yields an additional level structure, i.e. the 
one shown in Fig. [2] For clarity, we only show the dot levels that 
are located in the gap of the barrier regions and are, therefore, bound 
states. In the figure, we choose to present a situation with three bound 
states in total: Two of them are of the n = series (straight red lines 
in the center region) and a single one is of the n = - 1 series (straight 
green line in the center region), (b) Energy bands for a double dot 
setup. A single bound state (straight black line) is shown in the con- 
duction band of the left dot and two bound states are shown in the 
conduction band of the right dot. They are coupled via the contin- 
uum in the valence band of the central barrier which is enabled by 
the Klein paradox. 



wave functions of the left and right dots and e is the single- 
particle ground state energy. Note that the overlap integral 
vanishes if the states on the left and right dot belong to dif- 
ferent transverse quantum numbers q„ L + q llR . For the ground 
state mode, we have n^ — n« = 0, and the hopping matrix 
element can be estimated for d > L as 

t « 4sa 5* Wdzo,k exp(-d\k\), (8) 

where ao and 5o are wave function amplitudes (with dimen- 
sion 1/length) that are specified in App. C. As expected, the 
exchange coupling decreases exponentially with the barrier 
thickness, the exponent given by the "forbidden" momen- 
tum k in the barrier, defined in Eq. (0). The amplitude t 
can be maximized by tuning to a bound-state solution, where 
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FIG. 4: Ground-state wave function. Normalized squared wave 
functions \<¥ A \ 2 = \¥f\ 2 = \¥f } \ 2 and \¥ B \ 2 = \mf>\ 2 = \¥f } \ 2 for 
the bound state solution for the parameter choice e(Vbamer - V sate ) = 
0.6hvqo, and qoL = (n/3)L/W = 2 (indicated by the arrow in Fig.[2}. 
The corresponding energy is given by s « 1.31 • %vq . The dotted 
lines indicate the dot region < y < L. 



\s - eVbamerl approaches ftvqo (from below). Then, d\k\ < 1. 
Such a fine-tuning can be easily achieved in graphene quan- 
tum dots, where the small band gap allows to sweep through it 
and, therefore, to use conduction and valence band states (of 
the barrier region) to couple quantum dots. In Fig. Ob), we 
sketch the energy bands for the double dot case which shows 
how confined states in the two dots can be coupled via co- 
tunnelling processes through the continuum of states in the 
valence band of the central barrier region. Remarkably, this 
opens up the possibility for long distance coupling of electron 
spins because, in the limit |e - eVbamerl — > hvqo, the coupling 
t depends only weakly on the distance d between the quantum 
dots. However, already for bound state solutions with \k\d > 1, 
a coupling over a length exceeding several times the dot size 
is possible. For the situation where we couple two ground 
states in the quantum dots, we find, for instance, a solution, 
where \k\d = 4, d = 10 L, and the coupling can still be as 
large as t * 0.03e for highly localized qubits. This example 
is shown in Fig. [5] If we couple a ground state in the one dot 
with an excited state in the other dot, the hopping matrix ele- 
ment t can be even larger. The corresponding wave functions 
for that case are illustrated in Fig. [9] The values of t, U, and 
/ can be estimated as follows. The tunnelling matrix element 
t is a fraction of s as 30 meV (for a width of W ~ 30 nm), 
we obtain that t * 0.5 . . . 2.5 meV. The value for U depends 
on screening which we can assume to be relatively weak in 
graphene 12411 . thus, we estimate, e.g., U x lOmeV, and ob- 
tain J * 0.1 ... 1.5 me V. 

For the situation with more than two dots in a line, it turns 
out that we can couple any two of them with the others being 
decoupled by detuning. We mention here that our model is 
based on a single particle picture. Such a model effectively 
captures effects of Coulomb interactions as far as they can be 
described within the constant interaction model for quantum 
dots lT30tl . The Coulomb interaction then only shifts the energy 
levels in each dot by a constant. In Fig.|6j we illustrate the sit- 
uation of three dots in a line where the left and the right dot are 




d=10L 



FIG. 5: Long-distance coupling of two qubit ground states. The 

normalized squared wave functions |<F| 2 = \¥f \ 2 + \¥f >| 2 + \¥®\ 2 + 
lYg '| 2 of two qubits separated by a distance d = 10L, where L is the 
length of each quantum dot, are plotted next to each other. A ground 
state (of the series with the transverse quantum number n = 0) in the 
left dot is coupled to a ground state (of the same series with n = 0) 
in the right dot. The coupling is as large as t - 0.03e, where e is the 
ground-state energy. Furthermore, the qubits are highly localized, 
which can be seen from the ratio P(in)l P(out). Here, P{in) is the 
probability of the electron to be inside the corresponding dot and 
P(out) is probability to be outside the dot in the barrier regions. The 
parameters chosen for the potential (in units of Tivq a / e) are Vbameri = 

Vbarrier3 = 1, Vbarrier2 = 1.65, V gat el = V gat e2 = -0.5. 




FIG. 6: Triple quantum dot setup. The energy bands of a triple 
quantum dot setup are shown in which dot 1 and dot 3 are strongly 
coupled via cotunnelling processes through the valence bands of bar- 
rier 2, barrier 3, and dot 2. The center dot 2 is decoupled by detun- 
ing. The energy levels are chosen such that Ae2 <K Aei . The triple 
dot example illustrates that in a line of quantum dots, it is possible to 
strongly couple any two of them and decouple the others by detun- 
ing. This is a unique feature of graphene and cannot be achieved in 
semiconductors such as GaAs that have a much larger gap. 



strongly coupled and the center dot is decoupled by detuning. 
The tunnel coupling of dot 1 and dot 3 is then achieved via 
Klein tunneling through the valence band of the two central 
barriers and the valence band of the center dot. It is important 
for the long-distance coupling that the exchange coupling of 
qubit 1 and qubit 3 is primarily achieved via the valence band 
and not via the qubit level of the center dot - leaving the qubit 
state of dot 2 unchanged. Using the standard transition matrix 
approach, we can compare the transition rate of coupling dot 
1 and dot 3 via the continuum of states in the valence band of 
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the center dot (which we call Fvb) with the transition rate via 
the detuned qubit level of the center dot (which we call Fq B ). 
We obtain for the ratio (see App. D for the derivation) 



r V B/r QB ^(L/W)ln(4A/£ gap ), 



(9) 



where A * 6eV is the band width of graphene. There- 
fore, by increasing the aspect ratio L/W, it is possible to in- 
crease the rate Fvb with respect to Fqb- For L/W = 2 and 
£gap = 60meV, we find that Fvb /Fqb ~ 12, meaning that 
the qubit level in dot 2 is barely used to couple dot 1 and dot 
3. This is a unique feature of graphene quantum dots due to 
the small and highly symmetric band gap, which is not known 
to exist for other semiconducting materials. The availability 
of non-local interactions is important in the context of quan- 
tum error correction, as it raises the error threshold for fault- 
tolerant quantum computation [31]. In conclusion, we have 
proposed a setup to form spin qubits in quantum dots based on 
graphene nanoribbons with semiconducting armchair bound- 
aries. For such a system, we have calculated bound states of a 
tunable dot and outlined how two-qubit gates can be realized. 
We expect very long coherence times for such spin qubits 
since spin-orbit coupling and hyperfine interaction are known 
to be weak in carbon, see App. E. Furthermore, we have found 
that the high flexibility in tuning graphene quantum dots in 
combination with conduction band to valence band tunnelling 
based on the Klein paradox allows for long distance coupling 
of electron spins. Therefore, we propose a system which can 
serve as the fundamental building block for scalable and fault- 
tolerant quantum computing. 



APPENDIX A: GENERAL MODEL 

We now present in detail how to derive solutions for bound 
states in a graphene quantum dot. The dot is assumed to 
be rectangular with width W and length L as illustrated in 
Fig. 1. The basic idea of forming the dot is to take a strip 
of graphene with semiconducting armchair boundary condi- 
tions in jc-direction and to electrically confine particles in y- 
direction. Transport properties of a similar system have been 
discussed in Ref. B25II . 

The low energy properties of electrons with energy s in 
such a setup are described by the 4x4 Dirac equation 

-mvh\^ _^° + W(y)* = rfP. (Al) 



with the electric gate potential 



V(y) = 



'gate •> 
'barrier: 



(0 < y < L), 
otherwise. 



(A2) 



In Eq. ( lAU . cr x and o~ y are Pauli matrices , h is Planck's con- 
stant devided by 2n, e is the electron charge, and v is the 
Fermi velocity. The four component envelope wave function 
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t , i ; , , ' i - ' /( ) varies on scales large com- 
pared to the lattice spacing. Here, A and B refer to the two 
sublattices in the two-dimensional honeycomb lattice of car- 
bon atoms, whereas K and K' refer to the vectors in reciprocal 



space corresponding to the two valleys in the bandstructure of 
graphene. 

Plane wave solutions to Eq. (HJ take the form l29tl 
with 



X„l(x) = a n . 
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(A4) 



and 
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The complex number z n ,k is given by 

q„ + ik 



Z n ,k = ± 



(A5) 



(A6) 



The energy of the state in the barrier regions (y < and y > L, 
where V = Vbanier) is given by 



2V ba 



± hv 



(A7) 



In the dot (0 < y < L, where V = V gRte ) the wave vector k is 
replaced by k, satisfying 



s = eVgate ± hv yjql + k 2 . 



(A8) 



The ± sign in Eqs. dA6t - dA8b refers to conduction and va- 
lence bands. In the following, we concentrate on conduction 
band solutions of the problem (keeping in mind that there is 
always a particle-hole conjugated partner solution). 

The transverse wave vector q„ as well as the coefficients 
a n,±, a'„ ±, b„ y± , b' n ± of the n-th mode are determined (up to a 
normalization constant) by the boundary conditions at x = 
and x — W. We consider a class of boundary conditions for 
which the resulting parameters are independent of the longi- 
tudinal wave vectors k and k. We are particularly interested 
in semiconducting armchair boundary conditions defined by 
El 
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1 







-i2n/x/3 j 





x=W> 



(A9) 
(A10) 
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where yu = ±1 is defined by the width of the graphene strip 
W = «o(3M + jS), with M a positive integer («o = 0.246 nm is 
the graphene lattice constant), and 1 is the 2x2 unit matrix. A 
strip whose width is an integer multiple of three unit cells (yu 
0) is metallic and not suitable for spin qubit applications. The 
states of a semiconductor strip are non-degenerate (in valley 
space): 



w 



(All) 



with a„, ± = b' n ± = 0, a' n ± = b„ i± (for/z = -1) or a' n ± = 
b n ,± - 0, a„ i± - b' n± (for // = 1). Note that q„ determines the 
size of the gap for each mode n that is due to the boundary 
conditions. The size of the gap of mode n is given by 2hvq n . 
For concreteness, we consider the case of /u = 1 only. It can 
be shown that for the case of yU = -1, the bound states and 
normalized squared wave functions have exactly the same de- 
pendence on the parameters of a quantum dot, Vbamer, Vgate, 
W, and L. 

Our ansatz for a bound state solution at energy s to Eq. ([TJ 
then reads 



-iky 
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n,k 



if y < 0, 

if < y < L, (A12) 
if y > L. 



For bound states, the wave function should decay fory 
so we require that 



= iyjql-(.t 



■ eVbamer) 2 /(^) 2 



(A13) 



where fivq n > \s — eVbarrierl always has to hold. 

To find bound state solutions, we have to analyze the fol- 
lowing set of equations (coming from wave function matching 
at y = and y = L) 
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Zn,k 
1 

1 

Zn,k 



1 



1 



(A14) 



where a n = a'„a„-, f3„ = p' n a nA ,y „ = Y n a n _, S n = 6' n a n , + , and 
S = ikL. We can write Eq. (lA14t as 

-1 



Z n ,k 
1 



-z, 



~Z„i 



~z n - k 

-1 



\ 






Pn 




In 


/ 


fin j 



= 0. 



(A15) 



The allowed energy values s are readily determined by finding 
the roots of the determinant of the matrix on the lhs of the 
latter equation 



e 2S (Zn,k - Z„lY - (1 - Z n ,kZ„ lY = 0. 



(A16) 



A rather obvious solution is z„ % - 1 (which implies that 
S = 0) corresponding to 



s = ±hvq n + eVg ate . 



(A17) 



However, the corresponding wave functions to the allowed en- 
ergy solutions ( |A17t vanish identically. So, in order to pro- 
ceed, we have to find other (less trivial) solutions to the tran- 
scendental equation ( IA16b . Since the case z n ,k = Z„% only 
has the trivial solution ( IA17b . we can assume that z n ,k + Z n %, 
which means that Vbanier + Vgate- Then, we find that 



e s =± 



1 



' Zn,kZ n j; 



Zn,k 



(A18) 



^•n.k 



To analyze the solutions of Eq. ( 1A18I I. we distinguish two 
cases, one where 



S — ikL 



(A19) 



is purely imaginary (i.e., k is real), and another, where S is 
real. The two cases are distinguished by the criterion |e - 
eVgatel ^ hvq„ and |e - eV ea tel < fivq n , respectively. In the 
former case, since the lhs of Eq. jAl 8b has modulus unity, also 
the rhs must be unimodular, which is satisfied if in addition 
|e - <?Vbanierl - ^ wc ln- This case, where the equation for the 
argument of Eq. flA18t remains to be solved, is discussed in 
Sec.|B] The latter case where S is real has no solutions. 
Indeed, let us rewrite Eq. ( IA1 8b as follows: 



„2S 



1 



-2kk 



ql~(e- eVban-ierXe - eV g ate)/(^v) 2 + kk 



(A20) 



Taking into account that q\ > (s - eyban-ier)(s - eV gate )/(tivY 
and kk e R, we find that the left and right sides of this equation 
have different signs, therefore, Eq. ( IA18I ) has no roots for any 
purely imaginary k. 



APPENDIX B: BOUND STATE SOLUTIONS 

We now restrict ourselves to the energy window 

|e - eVgatel > hvq„ >\e- eVbamerl- (Bl) 

Then k is real, therefore |e s | = 1 and |z„j| = 1. Furthermore, 
z n ,k is real. We define z„ \ = e' e % where 



6„ = arctan(fc/q , „). 
It is easy to verify that in the energy window ( IB II ) 

1 — Zn,kZ n l 



= 1. 



Zn,k Z n k 

We can now rewrite Eq. (IA 18b as 

sin0„(l-z 2 ,) 



tan(feL) = 



2z„, k -(l + z 2 nl ) cos 6 n 



and further simplify this expression by using that 

k/q„ 



sin ft, = 



cos ft, 



yfl + (~k/q n ) 2 
1 

Jl + Ck/qnY 



(B2) 
(B3) 

(B4) 

(B5) 
(B6) 
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After some algebra, we obtain 



tan(JtL) 



-ikk 



(s - eVbamei-)(e - eV gate )/(hv) 2 - q 2 



(B7) 



The latter equation in combination with Eq. dAl 3b yields 
Eq. (7). Numerical solutions to Eq. ( IB7I ) are shown in Fig. 2. 

By applying different voltages to the gate and the barri- 
ers we shift the energy bands of the graphene ribbon un- 
der the barriers with respect to that of the quantum dot. A 
bound state in the quantum dot is allowed once the energy 
of the state hits the band gap of the barriers. If the differ- 
ence of the barrier and gate voltages AV = |Vban-ier - V ga tel 
is less than the energy of the gap 2hv\q n \ for n-th subband, 
k of a bound state lies in the interval [-k max , k max \, where 
£max is found from the condition s(k m - dx ) - Vbamer + hv\q„\ 
and, therefore, £ max = eAV -^/l + 2hv\q„\/eAV. The num- 
ber of bound states (for a given subband index n) is propor- 
tional to the length of the quantum dot L and is given by 
N = [kmaxL/n] (Jx\ is the integer just larger than x) . The 
number of the bound states of the n-th subband is maximal, 
when the barrier-gate voltage difference equals the energy 
band gap (AV = 2hv\q„\), and so 7V max = [ ^S\q n \L/n}. In 
the case of AV > 2hv\q„\, the top of the valence band of the 
graphene ribbon under the barriers becomes higher than the 
bottom of the conduction band of the quantum dot, therefore, 
there are no bound states with energies s < eVban-ier - hv\q„\ 
and \k\ of a bound state lies in the interval [fc m j n , £ max ], where 
£min is found from the condition s(k m i n ) = Vbanier - hv\q„\ 
(k m in = eAV - 2fiv\q n \/eAV), therefore, bound states lie in 
the energy window eVbamei- - M<7„| < s < eVban-ier + «v|<7nl (as 
shown in Fig. 2) and the number of the bound states is given 
by TV = {k max L/7T'] - {k min L/nl 

With increasing the barrier-gate voltage difference, 
a ra-th bound state appears at AVo = -hv\q n \ + 
%v -Jql + (n/L) 2 (m - l) 2 with the energy = V gat e + AVo + 
hv\q n \ and ends up at AV\ - hv\q n \ + fiv ^jq 2 + (nlL) 2 m 2 with 
the energy s$ = V gafe + AVi - %v\q n \ (see Fig. 2). 



APPENDIX C: WAVE FUNCTION 



0.06 



0.04 



0.02 




FIG. 7: Ground-state wave function. Normalized squared wave 

yTO|2 _ |VT/(«')|2 „„A IW_|2 



function |T A | 2 = IT™! 2 = |T\ A f and |T S | 2 



W>T\ 2 = \< 'I 2 for 



the ground state solution of a dot of length q„L = 5 with correspond- 
ing energy e = 1.101 hvq„. Here, e(Vbanier - Vgate) = 0.5 hvq n . The 
dotted lines indicate the dot region < yq„ < 5. 



(for fx — A,B each) finally determines the normalization con- 
stant of the wave function. With the ansatz dA12b we explicitly 
obtain the following components of the wave function 



¥f\x,y) 



-¥f \x,y) = 



a n z n , k e iq ' x e- ik y, 



iq n x -iky 



Pne knX e iky +7nZn . 
6 n e iq " x e ik(y - L \ 
a n e k '» x e- ik ?, 

p n z n ~ k e iq " x e ily +y n e i i» x e- ik * 
, 5 n Zn, k e iq " x e ik{ >- L \ 



iq n Xg-iky 



-iq„x„-iky 



Pne-^Jb +y„z n z< 

6ne -iq n x e ik(y-L) ^ 

' a n e- ic i" x e- ik y, 

p n z n ie-^ x e^ + y n e- iq » x e- ik *, 



(C3) 



(C4) 



(C5) 



(C6) 



-iq„x e ik(y-L) 



Following Brey and Fertig 112811 . we write the wave function 



y(x,y) 



-y ( P(x,y) 
-y ( P(x,y). 



(CD 



and give solutions for each component separately. As men- 
tioned above, the subscripts A and B refer to the two sub- 
lattices in the two-dimensional honeycomb lattice of carbon 
atoms and the superscripts K and K' refer to the two valleys 
in graphene. Note that the normalization condition rf28tl 

J ddx dy [pp«(jc y)\ 2 + WfHx, y)\ 2 ] = \ (C2) 



In the latter equations, the first line corresponds to the region 
in space, where y < 0, the second line to < y < L, and the 
third line to y > L. Thus, we obtain that 



\v;>(x, y r = pf7W)i 2 

\^f\x,y)\ 2 = ll-fW)! 2 . 



(C7) 
(C8) 



We now plot the normalized squared wave function of a 
ground-state solution and a excited-state solution of a dot with 
length q„L = 5 in Figs. [7] and [8] respectively. These are ob- 
tained from Fig. 2 under the choice that e( Vbanier - V gate ) = 
0.5hvq„. Evidently, the ground-state solution has no nodes in 
the dot region, whereas the excited-state solution has nodes. 
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FIG. 8: Excited-state wave function. 



function l^l 2 = l^l 2 



and \V B \ 2 = \¥ B K) \ 2 



Normalized squared wave 
\¥f >\ 2 for 

the first excited state solution of a dot of length q„L = 5 with corre- 
sponding energy e = 1.34 Tivq„. Here, e(Vbairiei - V g *te) = 0.5 Tivq,,. 
The dotted lines indicate the dot region < yq„ < 5. 





P(in)/P(out)=7.2 P(in)/P(out^=2.7 




^5arrier2 




1 ' ' 












l/gate 1 t=0.06£ v 

d=10L 





i 



is achieved via a continuum of states in the valence band of 
the barrier region as shown in Fig. 3(b) of the Letter. There- 
fore, the long-range coupling is enabled by the Klein paradox. 
In the weak tunneling regime, the hopping matrix element is 
given by 



(x, yy¥R(x, y)dx dy, 



(Dl) 



where *pL«(x,;y) = ^(x,y ± (d + L)/2) are the spinor wave 
functions of the left and right dots and s is the single-particle 
energy of the coupled levels. The integration in transverse re- 
direction is trivial and just gives a factor W. The integration 
in longitudinal y-direction can be restricted to the integration 
window y e [-d/2,d/2] if the wave functions are predomi- 
nantly localized in the dot regions. Then, the hopping matrix 
element can be estimated for d > L as 



t ~ 4easo6QWdzo,k exp(-d\k\), 



(D2) 



where ao and 6q are wave function amplitudes specified in 
Eqs. ( IC3b - dC6b . In Eq. dD2l) . we assumed that only levels 
of the series corresponding to the n = transverse mode are 
coupled. It is easy to relax this assumption because, if higher 
transverse modes form bound states, then only modes with 
«£ = Hr contribute to t, where is the transverse quantum 
number in the left/right dot. In Fig. [9] we demonstrate that 
a rather large coupling of t — 0.06e can be achieved over a 
distance as large as ten times the size of the quantum dots 
(see also Fig. 5 for comparison). Note that the qubits in this 
example are well localized in the corresponding dot regions: 
The probability of the electron in the left dot to be in the dot 
region P(iri) is 7.2 times larger than to be in the barrier regions 
P(out). For the right dot, the ratio of P(iri)/P(ouf) = 2.7 is a 
bit smaller but the electron is still predominantly localized in 
the dot region. 



FIG. 9: Long-distance coupling of a ground state and an excited 
state. The normalized squared wave functions pP| 2 = \¥f\ 2 + 
\¥f } \ 2 + \¥f>\ 2 + \¥f } \ 2 of two qubits separated by a distance 
d = 10L, where L is the length of each quantum dot, are plotted 
next to each other. A ground state (of the series with the transverse 
quantum number n = 0) in the left dot is coupled to an excited state 
(of the same series with n = 0) in the right dot. The coupling is 
as large as t = 0.06£, where e is the ground-state energy of the left 
dot. Furthermore, the qubits are still highly localized, which can be 
seen from the ratio P(in)/P(out). Here, P(in) is the probability of the 
electron to be inside the corresponding dot and P(out) is probability 
to be outside the dot in the barrier regions. The parameters cho- 
sen for the potential (in units of %vq /e) are V barricrl = V b arrier3 = 1> 
Vbamer2 = 1-65, V galcl = -0.5, and V gatc 2 = "0.9. 



APPENDIX D: LONG-DISTANCE COUPLING 

1. Long-distance coupling of two qubits 

Here, we discuss a particular example of long-distance cou- 
pling of two qubits separated by a distance d. The coupling 



2. Long-distance coupling in multiple quantum dot setup 

In Fig. 6, we propose a triple quantum dot setup in which 
dot 1 and dot 3 are strongly coupled and the center dot 2 is 
decoupled by detuning. It is important that dot 1 and dot 3 
are coupled via the valence band states of dot 2 and not via 
the (detuned) qubit level of dot 2. Otherwise, the spin of the 
decoupled qubit level would be affected by the coupling of the 
other qubits which is unwanted in the proposed long-distance 
coupling scheme. We assume that the gates that put the three 
dots in the Coulomb blockade regime are set in such a way that 
cotunnelling processes from dot 1 via dot 2 to dot 3 happen in 
the following order: First, an electron tunnels from dot 2 to 
dot 3 and then an electron tunnels from dot 1 to dot 2. The 
system is described by the Hamiltonian 

H = H + H T , (D3) 

where the kinetic term describes three qubit levels (a - 1,2,3) 
and the continuum of states in the valence band of dot 2 

H = ^ 2 Ea,o-al,o- a a,cr + ^ ^ 3fc,<r*£ .*lt > tr (D4) 
<*=l-3o-=m k (T=t,l 
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and the tunnelling Hamiltonian reads 



where 



H T = t ^ (a\ + a ] 2cr ai,cr + H.c.) 

cr 

+ 1 Yj i a ho- bk ^ + b l,cr a ^ + HX -)' (D5) 

In Eq. (ID4l l. a^o- and l>i jir annihilate electrons with spin cr 
in the qubit level of dot a and in the valence band of dot 2, 
respectively. We assume that E\ , T = #3 i<r (i.e. qubit 1 and 
qubit 3 are on resonance) and Aei = £2,0- - ~ £gap = 
2hvqo (see Fig. 6). In Eq. SD5\ . we make the approximation 
that all tunnelling matrix elements t depend only very weakly 
on energy and are real. The transmission rate from an initial 
state \i) to a final state |/) can be calculated using Fermi's 
golden rule 



W fi = ^\(f\T( Si )\i}\ 2 6(s f - £i ) 



(D6) 



with the transition matrix given by (up to second order in Hj) 

1 



T(s) = H T + H T 



s + irj - Hq 



H t + . 



(D7) 



We can put 77 = in the latter equation because we are only 
interested in off-resonant cotunnelling processes. This means 
that AS2 (see Fig. 6) should be finite because we want to 
have well-localized qubit states. The corresponding matrix 
elements of the T-matrix (lD71 i may be written as 

{T{e)) k # = (H T ) Kk , + Y(H T ) hk —(H T )v# +... (D8) 



k' 



8 - S k „ 



Now, we want to calculate |r I3 | 2 = |<3|r(J? 1)tr )|l)r, where |1> 
and |3) are the ground states of qubit 1 and 3, respectively. The 
lowest non-vanishing contribution of 1 7"i3 1 2 is of cotunnelling 
type, i.e., of fourth order in t. It is possible to separate the 
different contributions to IT^I 2 into three terms, namely 



<QB),2 
13 I 



+ ir 



13 I ' 



13 



(D9) 



In the latter equation, 



17]? | 2 determines the transition rate 

via the qubit level of dot 2 (the unwanted process), |7" 13 I 2 
determines the transition rate via the continuum of states 
in the valence band of dot 2 (the wanted process), and 



2 Re [(T® B) Y 



T{ 3 ] is the interference term of the two paths. 



It is straightforward to derive that (for a given spin <x of qubit 
2) 



T® B) (E hcr ) = 



£1,0- - E 2 ,c 



(D10) 



This has to be compared with 



r (VB) 



(El*-) 



r 



)E m /2 E UlT - E 
-^;ln(4A/ £gap ), 



(Dll) 



v (E) 



tJe 2 - (£ gap /2) 2 



(D12) 



is the density of states of the mode n — with £g ap = 2frvqo. 
In Eq. (ID1 II) . we integrate over the whole band width of the 
valence band (bounded by A « 6 eV). The approximate result 
in the second line of Eq. ( ID1 U holds for the hierarchy of ener- 
gies A » Zigap » E\ t(T . (In a more general case, the integral in 
Eq. dDl lb can still be evaluated analytically but yields a less 
compact expression.) 

The contribution coming from |r[j B) | 2 is evidently the 
smallest term of the three terms on the rhs of Eq. dD9l ). If 
we want to compare the rate that does not affect the spin of 
the qubit level in dot 2 (which we call Fvb) with the largest 
of the rates that does affect the spin of the qubit level in dot 2 
(which we call Fob), we can estimate 



Tvb 
Fqb 



n(VB),2 
13 1 



|2Re[(r;? B) )*r 1 ( !' B) ] 



-ln(4A/£ gap ). 



(D13) 



The latter result is Eq. (9). It shows that by increasing the 
aspect ratio L/W we can increase the weight of the coupling 
via the valence band states of dot 2 (which is wanted) as com- 
pared to the weight of the coupling via the qubit state of dot 2 
(which is unwanted). 



APPENDIX E: DECOHERENCE 

Finally, we give some arguments and rough estimates for 
the spin decoherence in graphene. It is generally believed that 
spin-orbit effects are weaker in carbon than in GaAs due to the 
lower atomic weight. Therefore, the dominating mechanism 
for decoherence will be the hyperfine coupling to the nuclear 
spins that are present in the material. 

The coherence time given by the hyperfine coupling can be 
estimated as 111 211 t a y/N/cA~ l , where is the number of 
atoms in the dot (for a typical graphene dot, « 10 4 ), cN 
is the number of atoms per dot with a nuclear spin, and A 
the coupling constant of the hyperfine interaction. Since we 
mainly deal with n orbitals in graphene, the contact hyperfine 
interaction is strongly reduced and the hyperfine interaction is 
dominated by its dipolar part. From the best available calcu- 
lations for the dipolar hyperfine matrix elements ll32l[33Tl . one 
obtains 

A dip = ^^/iUo.38 J ueV, (El) 

where yu is the nuclear magneton of I3 C, hb is the Bohr mag- 
neton, and /io the vacuum dielectric constant. This estimated 
value for Adi p is smaller than AcaAs ~ 90 //e V by more than 
two orders of magnitude. The natural abundance of 13 C is 
about c w 1% which yields, with the values quoted above, a 
coherence time of approximately t w 10 /is, about a thousand 
times longer than in GaAs. Unlike in GaAs, this value can be 
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improved by isotopic purification. Reducing the C content preselection of the dots without any nuclear spin to be used as 
by a factor of about 100 already decreases the average num- qubits. 
ber of nuclear spins per dot to about one. This allows for a 
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